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ABSTRACT 

Sloshing cold fronts (CFs) arise from minor merger triggered gas sloshing. Their detailed structure 
depends on the properties of the intra-cluster medium (ICM): hydrodynamical simulations predict the 
CFs to be distorted by Kelvin-Helmholtz instabilities (KHIs), but aligned magnetic fields, viscosity, 
or thermal conduction can suppress the KHIs. Thus, observing the detailed structure of sloshing CFs 
can be used to constrain these ICM properties. Both smooth and distorted sloshing CFs have been 
observed, indicating that the KHI is suppressed in some clusters, but not in all. Consequently, we need 
to address at least some sloshing clusters individually before drawing general conclusions about the 
ICM properties. We present the first detailed attempt to constrain the ICM properties in a specific 
cluster from the structure of its sloshing CF. Proximity and brightness make the Virgo cluster an 
ideal target. We combine observations and Virgo-specific hydrodynamical sloshing simulations. Here 
we focus on a Spitzer-like temperature dependent viscosity as a mechanism to suppress the KHI, but 
discuss the alternative mechanisms in detail. We identify the CF at 90 kpc north and north-east of the 
Virgo center as the best location in the cluster to observe a possible KHI suppression. For viscosities 
> 10% of the Spitzer value KHIs at this CF are suppressed. We describe in detail the observable 
signatures at low and high viscosities, i.e. in the presence or absence of KHIs. We find indications for 
a low ICM viscosity in archival XMM-Newton data and demonstrate the detectability of the predicted 
features in deep Chandra observations. 



1. INTRODUCTION 

Viscous stresses and thermal conduction in the intra- 
cluster medium (ICM) of clusters of galaxies play a crit- 
ical role in shaping the thermal structure of the cluster 
gas, in dissipating and redistributing the energy released 
by mergers and AGN outbursts, and controlling the feed- 
back between the ICM and the AGN inside the Bondi 
radius. In general, the transport processes in the gas 
must play a key role in the formation of structure down 
to scales of galaxies and star formation. 

The magnitude of the effective transport coefficients 
however is one of the longstanding and still unresolved 
questions. In an unmagnetized plasma, under the condi- 
tions of the ICM, both thermal conductivity and viscos- 
ity are l arge due to the strong temperature dependence 

T 5 / 2 dSpitzer||l956l |Sarazin||1988| ). For example, typi- 
cal Reynolds numbers in the ICM are 
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for length scales of A, velocities J7, particle density nicM 
and ICM temperature Ticm- In this equation we have 
included the suppression factor f ^ to parametrize an un- 
known degree of suppression of the viscosity. For a full 
Spitzer viscosity we have — 1. However, even weak 
magnetic fields strongly suppress diffusion of heat and 



momentum perpendicular to the field, and tangled mag- 
netic fields may lead to an overall suppression of conduc- 
tion and viscosity. The degree of suppression depends 
on the magnetic field structure and could be moder- 
ate (f fl ^ 0-3) or significant with « 1 (|Narayan fc| 



Medvedev 2001). The magnitude of the effective trans- 
port coefficients is relevant for numerous processes in 
galaxy clusters, e.g., the morphology and evolution of 



thus 



AGN inflated cavities ( Reynolds et aLf2005 Kaiser et al. 
2005j), the dissipation of AGN induced sound wa ves and 
balanci ng radiative cooling in cool cores dFabian 
2005a| ), and gas stripping from galaxies ( jNulsen 
On scales smaller than the coherence length, dif- 



et al 

fusion of energy and momentum should be anisotropic 
preferentially along the field lines. This can in itself 
lead to interesting effects and new instabilities like the 
magneto-thermal instabil ity and the heat flux driven 
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buoyancy instability (e.g.,|Balbus fc Reynolds 2010[ 

[Ku^[201Tp<unz et al.pOl^ . TKeTeP 

e processes for different aspects of ICM 



evance of these processes for different aspect 
physics is a topic currently being investigated by differ- 
ent groups. 

There are many theoretical arguments regarding the 
magnitude of the effective ICM viscosity, but the obser- 
vational evidence in the cluster ICM is controversial. For 
example, the coherent morphology of AGN ca vities can 
be explained by a substantial ICM viscosity (Reynolds 



et al. 



2005 Fabian et al. 2003), but turbulence in clus- 



ter cores has been invoked as a mechanism to amplify 
tiny p rimordial magnet ic fields to the values observed 
today (|Ryu et al.||2011|), and to broaden the heavy ele- 
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merit abundance peaks in cluster centers (Rebusco et al. 



2005, 2006 Roediger et al. 2007|). Several observations 



have been made to search for turbulence in the ICM. X- 
ray studies of the ICM in the Coma cluster estimated 
the power spectrum of density fluctuations in the gas by 
measuring th e spatial variation of the surface brightnes s 



measuring tne spatial variation ol tne surlace brigntnes; 
fluctuations (Sch uecker et al.|2 004; Ch urazov et al |2 012) 
Turbulence in the ICM is one possible origin of thes( 



possible origin of these 
fluctuations. Additionally, attempts to measure the tur- 
bulent velocity broadening of emission lines in clusters 
and early-type galaxies have been made with the XMM- 
Newton RGS ( jSanders et al.|[2010[ ). Future high spec- 
tral resolution observations will allow us to constrain the 
velocity power spectrum of the ICM by me asuring the 



centroids a nd widths of X-ray emission lines ( jZhuravleva 
etH1[2QT2 ). 



There are some places in galaxy clusters where we may 
directly witness the impact of the effective viscosity, cold 
fronts being one of them. Sloshing cold fronts (CFs) are 
density and temperature disco ntinuities in the ICM of 



galaxy clusters (see review by Markevitch & Vikhlinin 
2007, also for how they are distinguished from merger 
CFs). They are observed as arc-like edges in surface 
brightness, temperature and metallicity maps, wrapped 
around the cluster cores in spiral-like or staggered arc- 
like patterns (e.g., Figs. [I] and [3] later in this paper). They 
are attributed to sloshing or oscillatory motions of the 
ICM inside the cluster potentials, triggered by a merger 
that was too wea k to destroy the cluster core ( |Marke-| 
vitch et al .l|2001| |Ascasibar fc Markevit^ | 2QQ6p^uch 
sloshing CFs are commonly seen in clusters ( Markevitch| 
et al. 2003 Ghizzardi et al.|2010 ) and have b een observed 



Markevitch 



in detail m severa l clusters and groups (s ee 
fc Vikhlinin]|2007l [Roediger et al.||2012b| for' lists of ex- 
amples, also Sect. 6.3). 



Sloshing CFs are accompanied by shear flows along 
them and thus should be subject to the Kelvin-Helmholtz 
instability (KHI), as confi rmed by high resolution hy- 
drodynamical sim ulations (IZuHone et al.|2010l iRoediger 



|et al.||2011[ |2012b| or Fig. [2] later in this paper). The 
growtn of KHI mod es above a critical wav elength is sup- 
pressed by gravity ( |Chandrasekhar|[l961[ ) . The presence 
of ~ 10 to 50 kpc KHI modes in the hydrodynamical 
simulations demonstrates that this critical wavelength 
is sufficiently large to allow the growth of KHIs at ob- 
servable length scales (the same can be derived analyti- 
cally). Furthermore, KHIs can be suppressed at all wave- 
lengths by sufficiently strong magnetic fields aligned with 
the fronts ([Ch andrasekhar 1961; Vikhl inin fc Markevitch 



2QQ2||Keshet et al.|2010[JZuHone et al.|201l[ ), or at small 



wavelengths by a sufficiently strong viscosity. Conse- 
quently, the presence or absence of KHIs at sloshing CFs 
can be used to constrain these ICM properties. Exam- 
ples of both smooth CFs and disto rted ones, resembling 
the KHI, are k nown (see Sect. [673] for examples). While 
ZuHone et aT] ( |2010[ |2011[ |2012fl have studied the global 



impact of viscosity, magnetic fields and thermal conduc- 
tion on the sloshing process, a detailed attempt to derive 
the ICM properties from the observed structure of slosh- 
ing CFs has not been performed yet. This is the aim of 
this and a series of forthcoming papers. 

At a distance of only 17 Mpc (1"=80 pc) the Virgo 
cluster is the nearest galaxy cluster and can thus be ob- 
served with high resolution. The XMM-Newton mosaic 




Fig. 1. XMM-Newton MOS1+2 image of the ICM of the Virgo 
cluster, background subtracted and exposure corrected (0.5-3.0 keV 
band). M87 lies at the brightness peak. The sloshing cold front 
to the north (N) of the nucleus is denoted by the yellow arrows. 
It clearly continues to the north-west (NW) and somewhat more 
diffuse to the north-east (NE). 

around M87 in the center of the Virgo cluster displays 
a sharp edge in surface brightness at 90 kpc north-west 
(NW) of the clust er core as shown i n Fig. |T] Using a sub- 
set of this data, |Simionescu et , alT ( |2010[ ) demonstrated 
that this feature is a sloshing CF. A aet ailed analysis re- 
vealed a secondary and tertiary CF at 30 kpc and 20 kpc 



to the south-east and NW, respectively. Roed iger et al 
( |2011[ ) simulated ICM sloshing specifically for the Virgo 
cluster and inferred the most likely merger scenario from 
a detailed comparison to the observations. The proximity 
of Virgo, the high X-ray surface brightness and favorable 
merger geometry in the plane of the sky allow us to probe 
the shear flow interface with a spatial resolution of ~100 
pc. 

We have re-simulated the gas sloshing in Virgo at high 
resolution, with different levels of ICM viscosity as de- 
scribed in Sect. [2] We demonstrate that even a viscosity 
of about 10% 01 the Spitzer value efficiently suppresses 
the KHI in Sect. [3J and compare the simulation results 
to an analytical estimate. Sect. [4] describes observable 
features at low and high viscosity. We report marginal 
evidence for a low ICM viscosity in t he e xisting shal- 
low XMM-Newton observation (Section 5.1) of the Virgo 
cluster CF. We demonstrate that the presence or ab- 
sence of KHIs can be detected with deep Chandra ob- 
servations. We discuss features which could be detected 
with very deep observat ions with either Chandra or fu- 
ture telescopes (Section 5.3|). We discuss model uncer- 



tainties, alternative KHI suppression mechanisms, and 
briefly compare to other CF clusters in Sect. [6] Section [7] 
summarizes our results. 

2. METHOD 

We perform idealized minor merger simulations where 
a gas free subcluster of 2 x 1O 13 M moved through the 
Virgo cluster from west to east, passing 100 kpc south 
of the cluster core 1.7 Gyr ago. |Roediger et al. (2011) 
inferred this merger configuration as the most likely one 
from a detailed comparison of observed and simulated 
CF properties including orientation, distance to the clus- 
ter center, and contrasts in brightness and temperature 
across the fronts. Our simulations use the FLASH code 
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(version 3.3, Dubey et al. 2009) and utilize the same 
method and initial conditions as described in Roediger 
let al.l (|2011|) unless stated otherwise. 



2.1. Resolution 

We resolve the cluster core with a grid cell size of 0.5 
kpc out to 110 kpc from the cluster center. The FLASH 
code captures the growth of the KHI even at a resolution 
of 16 grid cells per perturbation scale length (Roediger et 
al., in preparation, also Fig. 11 in the Appendix), hence 
we can resolve modes down to~8 kpc. We performed the 
same simulations at a coarser resolution of 1 kpc and 
found consistent results, but present the higher resolu- 
tion results here. Further snapshots of a resol ution test 
in the inviscid case can be found in Fig. Al of Roediger 
Convergence is to be expected for sim- 
ligh viscosity, because viscosity damps the 



et al.| ([2011). 

ulations at 

instability at small wavelengths first. Thus, as long as 
the smallest unstable wavelength is resolved, increasing 
the resolution does not bring any new insights. At zero 
viscosity the simulations cannot formally converge, be- 
cause the KHI grows fastest at small wavelengths. How- 
ever, in the presence of multi-wavelength perturbations 
the largest modes which had sufficient time to grow dom- 
inate the visual appearance of a KH unstable interface. 
The presence of smaller modes may add to the widening 
of the interface, which is captured in our simulations. 
Hence, our simulations sufficiently resolve the observable 
features. Naturally, the exact patterns of individual KH 
rolls will depend on the details of the seeding, which is 
inaccessible in any case. 

2.2. Viscosity 

We include a temperature-dependent dyna mic vis- 
cosity /i = 5500gcm-VV M ( T iCM/10 8 K) 5 / 2 ( ISarazin 
1988, assuming a Coulomb logarithm of 40), where we 
parametrize the unknown level of viscosity suppression 
by the suppression factor f fl . Here our simu l ations differ 
from the ones presented by |ZuHone et al.| ( |2010[ ), who 
used a constant viscosity throughout the cluster. This 
may overestimate how effectively viscosity suppresses the 
KHI because it neglects the lower viscosity on the cooler 
side of the CF that would apply for a Spitzer-type viscos- 
ity. As pointed out above, the nature of the effective ICM 
viscosity is still unkn own, and both models are valuable. 
In fact, the choice of |ZuHone et al. (2010) was optimal 
for the focus of their work, which was to demonstrate 
the inability of viscous gas sloshing to prevent cooling 
flows, even in the presence of the highest reasonable vis- 
cosity. Using a temperature dependent viscosity provides 
a more stringent test of its ability to suppress the KHI, 
since having a lower viscosity on the cool side of the CF 
reduces viscous suppression of the KHI. Thus, this model 
is the better choice for our current purpose. Moreover, 
th e Virgo cluster i s muc h cooler than the cluster studied 
by ZuHone et al. (2010), which might further reduce its 



effectiveness at suppressing KHI. 

A full Spitzer thermal conductivity would erase tem- 
perature discontinuities rapidly: at a temperature of 
3 keV and density of 3 x 10 -3 cm -3 , the conduction 
timescale for temperature gradients with scale lengths 
of 3 k pc, i.e. across the fronts, is about 0.3 Myr ( |Sarazin| 
|1988[ ), much shorter than the age of CFs. Even if ordered 



magnetic fi elds prohibit the cond uction directly across 
the fronts, ZuHone et al. (2012) have shown that full 
Spitzer conductivity still erases the temperature discon- 
tinuity at the CFs in a few hundred Myr. The presence 
of temperature discontinuities at the observed CFs im- 
plies strongly suppressed thermal conduction, hence we 
neglect it. 

2.3. KHI seeding 

The KHI needs seeds, initial perturbations, to grow. In 
contrast to earlier work we do not rely on seeding by the 
discretization of the computational grid, but introduce 
intentional seeds by repeatedly adding small perturba- 
tions in the north-south velocity every 100 Myr, from 
pericenter passage on. The perturbation velocity is a su- 
perposition of sinusoidal waves of wave lengths 40, 20, 10, 
. . . , 1.25 kpc with a total amplitude of lOkms -1 at the 
first perturbation event. The amplitude is constant with 
wavelength. The amplitude decreases linearly at sub- 
sequent perturbation events such that the perturbation 
is zero at 2 Gyr after pericenter passage. In this man- 
ner, we introduce a perturbation of a few percent or less 
of the typical shear velocity at each perturbation event, 
but keep the total perturbation velocity well below the 
shear velocity. Perturbations of this level or higher can 
arise e.g. from ICM tu rbulence due to AGN activity or 



the motion of galaxies ( |Roediger fc Bri iggen 2008 ; |Vazza 
|et al. 2012). While our perturbations are highly ideal- 
ized in their power spectrum and direction, they serve 
the purpose as seeds for the KHI. 

2.4. Simulation runs 

We simulate the merger history from 1 Gyr prior to 
pericenter passage to the fiducial time at 1.7 Gyr after 
the pericenter passage for the low viscosity case of / M = 
10 -3 . Here KHIs occur from 0.5 Gyr after the pericenter 
passage on. We re-simulate the final 0.5 Gyr, i.e. from 1.2 
Gyr after the pericenter passage on, at higher viscosity 
with = 0.01 and 0.1. Besides reducing computational 
expenses, this procedure poses a more stringent test on 
the ability of the viscosity to suppress KHIs because the 
instability is clearly developed at 1.2 Gyr. Additionally, 
we performed an inviscid control run. 

3. SUPPRESSION OF KHI BY VISCOSITY 



Junk et al. (2010) derived a dispersion relation for 
the viscous KHI by standard linear perturbation anal- 
ysis for the case of constant kinematic viscosity. It pre- 
dicts reduced growth of the KHI for Reynolds numbers 
< / (Pcoid+Phot) • e £ or Reynolds numbers smaller 

V PcoldPhot ' J 

than ~ 10 for modera te density contrasts of around 2 
( |Roediger et al.||2012a ). The shear velocity and the per- 
turbation length scale enter the Reynolds number as the 
relevant velocity and length scale (Eqn. [TJ. This agrees 
with the general physical expectation that viscous forces 
are equally important to inertial forces at a Reynolds 
number around 1. Thus for a ratio of 2 between the 
density on the cold side of the northern front, p C oid, and 
that on the hot side of the front, phot? m Virgo, a full 
Spitzer viscosity should affect perturbations of length 
scales smaller than ~ 10 kpc. 

Figure [2] demonstrates the effect of the viscosity in the 
simulations. It displays temperature slices in the orbital 
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plane at the final timestep for = 0, 10 3 , 0.01 and 0.1 
from top to bottom. At a viscosity < 10 -3 , all CFs 
are clearly distorted and made ragged by the KHI. With 
increasing viscosity, the fronts become less ragged, and 
structures at progressively larger scales are suppressed. 
Interestingly, even the small physical viscosity of 10 -3 
Spitzer erases some of the smallest perturbations present 
in the inviscid simulation. Finally, the fronts are almost 
completely smooth in the high viscosity case (/ M = 0.1) 
except for two large KH rolls separated by ~ 40 kpc along 
the SW, where the shear flow is strongest (~ 500 kms -1 ). 
Distortions at smaller length scales though are absent at 
high viscosity at this location as well, whereas smaller 
distortions are present at this location at lower viscosity. 



Our simulations demonstrate that the viscosity is more 
efficient in suppressing the KHI than expected from the 
linear analysis. There are several reasons for this dif- 
ference: the growth time is derived from the linear sta- 
bility analysis. The effect of the viscosity is to reduce 
shear velocities, which applies to the flow parallel to the 
interface as well as the velocity perturbation in the per- 
pendicular direction. Thus, while the linear analysis pre- 
dicts only a slowed growth of the KHI, but still a growth, 
at longer timescales viscosity should shut off the growth 
completely and thus be more efficient than expected. We 
demonstrate this behavior analytically and numerically 
in a separate publication (Roediger et al., in prepara- 
tion). Furthermore, the sloshing CFs are curved inter- 
faces embedded in a background gravitational potential 
and a stratified atmosphere, whereas the analytic esti- 
mate assumes a planar interface, no stratification and no 
gravity. The gravity at the northern CF suppresses KHIs 
at wavelengths above ~ 50 kpc and will thus slow down 
the growth of instabilities at somewhat smaller wave- 
lengths, i.e. the wavelengths seen in the simulations. Fi- 
nally, the sloshing CFs are special contact discontinuities 
that are continuously reformed by the sloshing process, 
which interacts with the KHI and may modify its growth 
at long timescales. For example, the outwards motion of 
the CF stretches the relevant waveleng ths, which reduces 
the g rowth rate and amplitude ( jChurazov 
2004| . 
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4. OBSERVABLE FEATURES 

4.1. X-ray images 

We calculate synthetic X-ray images by projecting 
2 icm^(^icm) along the line-of-sig ht (LOS), where A(T) 



is the cooling function according to [Sutherland fc Dopita 
(119931) and we assume a metallicity of 0.3 solar. Figure [3] 
displays the predicted X-ray images for the viscosity sup- 
pression factors fn = 10 ~ 3 and f ^ = 0.1. We will refer 
to these two viscosities as low and high viscosity, respec- 
tively. 

As the inner CFs are likely to be distorted by the AGN 
activity in the Virgo cluster center ( Forman et al.|20'07 ), 
we focus on the structure of the outer northern front. 
The shear flow along this front is weakest in the NW, 
leading to a smooth, sharp front in the NW independent 
of viscosity. Along the north (N) and the east of the front 
the shear flow is stronger (~ 300 kms -1 ) and forms dis- 
tinct structures depending on the viscosity (Fig. [3|. At 
high viscosity, the front forms a smooth arc here as well, 




od>| / A 

Fig. 2. — Temperature slices in the orbital plane at the final 
timestep, for Spitzer-type, i.e. temperature dependent, viscosities 
with suppression factors = 0, 10 _3 ,0.01 and 0.1 from top to 
bottom. Increasing the viscosity erases progressively larger sub- 
structure along the fronts. We have oriented the images such that 
they compare to the situation observed in Virgo, i.e. north is up 
and west is right (see Roediger et al. 201 1] for details). 



but it is ragged at low viscosity. Individual KH rolls 
at ~ 15 kpc size can be identified as triangular-shaped 
irregularities. They give the front a saw tooth like ap- 
pearance. A further characteristic pattern is multiple 
adjacent brightness edges with a spacing of about 5 kpc 
parallel to the main front. We have labelled these fea- 
tures in the zoom-in in the top panel in Fig. [3] 
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Fig. 3. — Synthetic X-ray images in arbitrary logarithmic scale 
for the low (/^ = 10 -3 ) and high viscosity (f^ = 0.1) case in the 
middle and bottom panel. The top panel is a zoom-in in linear 
brightness scale on the northern part of the outer CF as marked 
by the dashed rectangle. Prominent structures are labelled. 



4.2. Brightness profiles 

Even though the CFs are discontinuities in density and 
temperature, they do not lead to true discontinuities in 
surface brightness, but in its slope. The X-ray brightness 
profile of a spherical cloud with a power-law density pro- 
file can be approxima ted by Sx (d) oc >/— d near its edge 
dVikhlinin et al |2QQ1[ ). In our notation, the coordinate d 
increases with increasing distance from the cluster cen- 
ter and is defined to be zero at the CF, i.e., — d is the 
distance to the projected edge measured inwards. The 
minus-sign is required because the brightness decreases 
with increasing distance to the center. Figure [4] compares 
brightness profiles across the CF towards the north-east 
(NE) for low and high viscosity. We extract profiles in 
narrow, 5° wedges, corresponding to a width of 9 kpc 
along the front, and for a wide wedge of 60° (see top panel 
in Fig. H). The vertex of the wedges is not the cluster 
er, but 
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Fig. 4. — X-ray surface brightness profiles across the outer CF 
towards the NE. The extraction wedges are shown in the top panel. 
The individual wedges 1 to 12 are 5° wide. The first line plot com- 
pares profiles for the total 60° wedge in the high and low viscosity 
case, see legend. For clarity, we have offset both cases by 0.1. We 
fit each profile by the expected brightness profile for a broken power 
law (Eqn.[2j. Thick lines are for the simulation data, thin lines are 
the corresponding fits. The bottom three panels show the same 
comparison for selected narrow wedges as labelled in each panel. 



center, out we aligned the annuli of the wide wedge with 
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the CF. While this makes the detailed interpretation of 
the profiles more difficult in terms of assumed geometry, 
this choice optimizes the sharpness of the edge in the 
profiles, which is the main purpose here. In all cases the 
CF appears as a steeper slope in the brightness profiles. 
We fit the profiles within 15 kpc around the edge with 
the expected broken power-law profile. More precisely, 
we allow the position of the edge dcF to be a free pa- 
rameter by working in the offset coordinate x = d + dcF 
along the profile, and include a linear background profile 
outside the CF: 

S x (d) = S(x - d C F) 

= {A y/(-x + d C F)/kpc + So if x < d C F, / 2 ) 
\ m(—x + dcv) + So if x > d C F, 

The amplitude of the jump Aq, the position of the edge 
dcF, the background brightness So, and the slope of the 
background profile m are free parameters. We have per- 
formed the same exercise for profiles towards the NW 
(not shown here), where the CF appears undistorted in 
the surface brightness images, and find good agreement 
with the expected >/— d shape. In the high viscosity case, 
profiles across the NE edge agree well with this expecta- 
tion as well (Fig. [4]). In the low- viscosity case, however, 
profiles across the ragged NE edge show a distinct de- 
viation from the >/— d profile. The slope at the edge is 
shallower, the edge appears washed or blurred out. In 
many of the narrow wedges, the profiles display multiple 
edges, corresponding to the multiple adjacent fronts seen 
in the image. 

5. CURRENT AND FUTURE OBSERVATIONS 
5.1. Evidence for KHIs in existing data 

The northern CF in the Virgo cluster, as evident in the 
XMM-Newton mosaic around M87, is highlighted by ar- 
rows in Fig. [T] While this CF appears sharpest, i.e. has 
the steepest brightness gradient, in the N-NW, it con- 
tinues in a somewhat diffuse mann er to the NE (see als o 
brightness residual map in Fig. 1 in Roedig er et al.|20lT ). 
This matches the predictions of our simulations at low 
viscosity. Due to the absence of a shear flow along the 
NW, the CF should be sharp here. Along the N and 
NE, our simulations show a stronger shear flow and con- 
sequently KHIs at low viscosity. In this shallow XMM- 
Newton observation, individual KH rolls cannot be de- 
tected, nor their absence confirmed, because of several 
factors. These include the broad- winged PSF of XMM- 
Newton, the limited counting statistics of the data (6.5 
ks of good time for the PN camera after flare filtering 
in this region of the CF), and the complex non-uniform 
detector response and background. 

To quantitatively confirm the visual impression of the 
smeared-out front in the NE we extracted surface bright- 
ness profiles from the archival XMM-Newton data across 
the CF both in the NW and the NE. Background was 
estimated using the dark sky background files, and the 
ODF files were reprocessed to incorporate the most up 
to date calibration data and filtered for periods of high 
background. We chose profile wedges which cover about 
5 arcmin (24 kpc) along the CF and carefully aligned the 
annuli of these wedges with the CF such that the radius 
of curvature of the sector matched that of the CF. The 
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Fig. 5. — Surface brightness profiles across the NW (green/light 
grey in print) and NE (black) regions of the Virgo cluster sloshing 
CF, extracted from the archival XMM-Newton data. The zero 
point on the x axis denotes the approximate position of the surface 
brightness edge of the CF. The scaling of the surface brightness is 
arbitrary. The profile wedges cover about 5 arcmin or 24 kpc along 
the CF. Note the sharp rise of the surface brightness behind the 
CF in the NW and the slower rise at in the NE. 

profiles are shown in Fig. |5j In the NW, the brightness 
profile exhibits a steep rise as expected for an underlying 
density discontinuity. In the NE, however, the bright- 
ness rises less sharply, indicating a smeared out front as 
expected for a low ICM viscosity. A similar apparent 
smearing may arise from a mis- alignment of the profile 
annuli with the front. We experimented with different 
extraction regions and could not get a sharper edge in 
the NE than shown in Fig. [5] Hence, the blurred CF in 
the NE suggests the presence of KHIs and thus strongly 
suppressed viscosity in the ICM. Individual KH rolls may 
be visible in a deep Chandra observation. 

Our simulated X-ray images exhibit two prominent KH 
rolls in the SW of the cluster core at both high and 
low viscosity. If present in this evolutionary stage, they 
might be detectable in the current XMM-Newton image. 
However, the exact shape and height of the KH rolls de- 
pends on the details of their seeding, which includes a 
full turbulent spectrum and perturbations from the cen- 
tral AGN outbursts in the real Virgo cluster, which is 
more complex than in our simulations. Furthermore, the 
rolls are less clear at a somewhat less optimal orientation 
of the LOS with respect to the merger plane. Thus, the 
basic prediction of the simulations is that the CF in the 
SW is especially susceptible to the KHI, and the XMM- 
Newton image indeed shows a less sharp front in this 
region. We note that these large KH rolls appear less 
clear in the X-ray images at low viscosity, because their 
borders are washed out by smaller-scale instabilities. At 
a viscosity of 0.1 Spitzer, these rolls appear much clearer 
and might be visible in the shallow XMM-Newton data, 
but they are not. Instead, the CF in the SW simply ap- 
pears blurred out over ~ 20 kpc, as would be expected 
at low Virgo ICM viscosity. Thus, the structure of both 
the NE and the SW fronts suggests a strongly suppressed 
ICM viscosity. 

5.2. Detectability of Kelvin- Helmholtz Instabilities in 
deep Chandra observations 

5.2.1. Direct imaging 
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Fig. 6. — Simulated X-ray images of the northern sloshing CF in the Virgo cluster at different viscosities. The top and bottom rows 
are for low and high viscosity (10 -3 and 0.1 of the Spitzer value), respectively. The left-hand-side column shows noiseless images, in 
the right-hand-side column we added a random Poisson deviate to match the surface brightness and noise level of a simulated 300 ks 
Chandra/ ACIS-I observation. The structure of the CF differs between low and high viscosity. The KHIs can be clearly seen in the former 
case (see labels), in both the ideal and in the noisy image. 



The left panels in Fig. [6] show the direct comparison of 
surface brightness images at high and low viscosity along 
the northern sloshing CF of the simulated Virgo cluster, 
i.e. a smooth front in the high viscosity case and a ragged 
front at low viscosity. This field of view corresponds to 
two ACIS-I pointings of the Chandra X-ray observatory. 

To evaluate the detectability of these structures in real 
observations, we match the count density in the simu- 
lated images to the surface brightness measured for the 
Virgo cluster CF in the XMM-Newton exposure. Using 
PIMMS, we scale it to a 300 ks Chandra /ACIS-I observa- 
tion. With this exposure time, there will be ~100 source 
counts per 0.5 kpc x 0.5 kpc pixel (6"x6") just behind 
the CF. A random Poisson deviate is then added to each 
0.5 kpc x 0.5 kpc pixel to simulate the noise of a real ob- 
servation. The resulting noisy images are shown in the 
right column of Fig. [6] We neglect background in these 
idealized simulations, because the count rate from the 
background is <10% of the rate from the gas inside the 
cold front and thus insignificant. 

The KH rolls in the low viscosity case are clearly visi- 
ble in the data with the random deviate added, and they 
are again absent in the high viscosity case. The KH rolls 
span spatial scales of several kpc. The variations in sur- 
face brightness due to the KH rolls in the low viscosity 
simulation are as large as ~20%. There are ~4000 counts 
in a 5 kpc x 2 kpc region behind the cold front in the 
simulated Chandra image. Thus we could in principle 
detect variations in surface brightness of ~5% at 3a con- 
fidence with a real observation. The KH rolls on scales of 
a few kpc, if present, would be easily detectable at more 
than 10a confidence in such a deep Chandra observation. 
Additionally, we could observe KH rolls on smaller scales 
at lower significance. Twenty percent variations of sur- 
face brightness could be detected at 4a confidence in 1 
kpcxl kpc scale regions. 



5.2.2. Profiles 

To further demonstrate the detectability of the KH 
rolls in the simulated data, we extracted surface bright- 
ness profiles across the CF in 1.5° wedges in both data 
sets with the random Poisson deviate added. Two ex- 
amples are shown in Fig. [7| We follow the classic ob- 
servational data analysis strategy and fix the vertices of 
the wedges at the cluster center. The wedge opening an- 
gle of 1.5° corresponds to a linear distance of 2.5 kpc or 
0.5 arcmin at the CF and is thus much narrower than 
in our analysis of the XMM-Newton data. Despite the 
added Poisson noise the predicted multiple edges in the 
low viscosity case can be clearly detected and are marked 
by vertical lines. The spacing between the surface bright- 
ness edges (i.e. the spacing between the vertical lines - 

3 to 4 kpc) is roughly the height of the KH rolls. The 
height of these rolls is typically a third or half of the scale 
length of the KH rolls. Thus the edges spaced roughly 

4 kpc apart in the low viscosity profile signify KHI scale 
lengths of ~10 kpc. The presence of these edges in a 
real observation would immediately give the typical scale 
length of KH rolls, and these rolls should be present if 
^icm <^ ^spitzer- In contrast, the surface brightness pro- 
file of the high viscosity simulation can be well described 
with a single power law with no evidence of an edge or 
change in slope inside the contact discontinuity of the 
CF. 

5.2.3. Statistical characterization 

At low viscosity, there is a wealth of structure inside 
the CF. We created the power spectral density (PSD) 
function of the emission behind the CF in both cases to 
search for a signature in the Fourier domain of the spa- 
tial scales where the KH rolls are present. The PSD is 
indeed larger in the low viscosity simulation on scales 
where the KH rolls are present, but we could not di- 
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Fig. 7. — Surface brightness profiles in 1.5° wedges (0.5 arcmin 
or 2.5 kpc at the CF) across the cold front in the low (top - 10 -3 
Spitzer) and high (bottom - 0.1 Spitzer) viscosity cases. The CF 
is at ~90 kpc. Note that the surface brightness profile of the low 
viscosity simulation shows multiple edges denoted by the vertical 
dashed lines. Such discontinuities are indicative of density jumps. 
They arise due to the KH rolls. The surface brightness profile for 
the high viscosity case displays only a single steep edge, indicating 
the absence of KH rolls. 



rectly recover the spatial scales at which the KHIs were 
induced. We can confirm the presence of the KH rolls 
using the PSD, but we do not expect to be able to ex- 
tract the specific spatial scales of the KH rolls in Fourier 
space. As an alternative, we produced brightness pro- 
files along the CF and derived their ID power spectra. 
In a narrow annulus along the CF the KH rolls appear 
as ripples in the azimuthal profile. To collect a sufficient 
number of counts, the annulus needs a certain width, 
but increasing the width of the annulus tends to aver- 
age out the surface brightness wriggles due to the saw- 
tooth pattern of the KH rolls. As a consequence, the 
ID power spectrum analysis does not reveal more infor- 
mation than can be obtained by visual inspection of the 
images. The Fourier space analysis is generally better 
suited for (quasi-) periodic signals, but we have only a 
few KH rolls in the simulated ACIS-I image. 

5.3. Structures in very deep observations 



The differences between high and low viscosity sloshing 
are more clearly demonstrated using an edge detection 
algorithm. We applied a Sobel filter, a discrete differen- 
tiation operator often used for edge detection, to both 
simulated surface brightness images as shown in Fig. [8j 
These filtered images were created using the raw (i.e., 
noiseless) simulated image and represent very deep ob- 
servations by Chandra or a future large area, low back- 
ground X-ray observatory where the noise from counting 
statistics is small. There are several features of note. 
First, the sloshing cold front is a sharp, continuous arc 
in the high viscosity simulation, whereas the CF is dis- 
torted in the low viscosity case. The approximate scales 
of the KHIs can be clearly seen in the filtered, low vis- 
cosity image. Second, the flow of the gas interior to the 
CF in the low viscosity case is turbulent, evident from 
the wealth of structure inside the cold front. Even in 
the presence of significant Poisson noise, these surface 
brightness fluctuations may be detectable in the power 
spectral density of fluctuations of the surface brightness 
in this region. 

6. DISCUSSION 

6.1. Model uncertainties 

In this section, the impacts of various assumptions on 
the conclusions are discussed. 

6.1.1. Merger history 

The cluster model and the assumed merger history 
determine the dynamical parameters relevant for the 
growth of the KHI at the CFs, namely the local density, 
temperature, their contrasts across the front, and the 
amplitude of the shear velocity. The density, tempera- 
ture and their contrasts across the CF are constrained 
observationally and matched in our simulation, hence 
these quantities are sufficiently accurate. The amplitude 
of shear velocity will depend somewhat on the merger 
strength. However, our results regarding the growth or 
suppression of the KHI are not sensitive to the exact 
shear velocity within a factor of 1.5. The shear velocity 
is unlikely to be outside this interval, because this im- 
plies either a much stronger or much weaker merger. The 
former would lead to larger contrasts of brightness and 
temperature across the fronts than observed and an over- 
all front distortion of the sloshing structure. The latter 
would lead to too weak contrasts or even the absence of a 
true discontinuity at the northern front ( |Roediger et al.| 
2011). Neither case is observed. Moreover, a viscosity of 
only 0.1 of the Spitzer value suppressed the KHI along 
the northern CF in our fiducial merger. A full Spitzer 
viscosity will certainly suppress the KHI in a moderately 
stronger merger. Hence our conclusions are little affected 
by uncertainties in the merger history. 

6.1.2. Line-of-sight (LOS) 

For our analysis we have assumed a favorable LOS, 
i.e. perpendicular to the orbital plane. The spiral-like 
appearance of the Virgo CFs favors a LOS close to this 
configuration. To test the dependence of the reported 
features on the merger geometry we have produced syn- 
thetic X-ray images for different LOSs. To enhance the 
relevant features, we have converted the X-ray images 
into brightness residual maps by dividing each image by 
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Fig. 8. — Simulated images of the Virgo cluster CF with Sobel filter applied, low viscosity on the left, high viscosity on the right. This 
filter selects regions of steep brightness gradients, i.e. edges. Note the difference in the structure along the CFs in the two images, and the 
clear evidence of turbulent flow behind the CF in the low viscosity (i.e., high Reynolds number) case. 



its azimuthal average, 
shown in Figs. [9] and 
xly 



The resulting residual maps are 
The high and low viscosity 
LOSs. At high viscosity, the 
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case differ clearly for 
outer CF always appears as a smooth arc with a sharp 
edge, whereas the edge is blurred out by the KHI at low 
viscosity, and the CF appears ragged and fuzzy. How- 
ever, if the merger is close to the plane of the sky, the 
front in the NW is sharp, even at low viscosity. The fact 
that the observed northern front in Virgo shows a rather 
sharp edge in the NW supports the conclusion that the 
merger in Virgo occurred almost in the plane of the sky, 
as suggested by the morphology of the brightness resid- 
ual maps. In this case, the best location to study the 
difference between existing and suppressed KHIs will be 
the front between N and NE, not in the NW. 

6.1.3. Resolution and seeding 

At zero viscosity and with ever increasing resolution, 
the simulations should show progressively smaller KHIs. 
Naively, one might expect the KHIs to fully smooth out 
the CF. However, the presence of KH rolls at a given 
timestep depends on the seeds that have been available 
as well. A full parameter study on this subject is be- 
yond the scope of this paper, but we stress that e.g. in 
clouds in the Earth's atmosphere, KH rolls of different 
wave lengths can be observed simultaneously without the 
cloud boundary being fully blurred out. In our simula- 
tions we simultaneously observe KH rolls of length scales 
differing by a factor of a few as well. Hence we conclude 
that individual KH rolls may be observable with suffi- 
cient spatial resolution. Moreover, a blurring out of the 
CF edge should be clearly detectable in the brightness 
profiles as shown above. 

6.1.4. Simplification of viscosity 

Our simulations describe the momentum diffusion in 
the ICM with an effective viscosity that includes the ex- 
pected Spitzer temperature dependence. Nature may be 
more complex than this. The effective ICM viscosity on 
scales of > 0.1 kpc arises from a complicated interplay of 
the plasma with the magnetic field. Various processes, 



e.g. the gas sloshing itself, impact the amplitude and ge- 
ometry of the magnetic fields in the cluster core, and the 
effective viscosity could vary spatially and temporally. 
Furthermore, due to the tiny electron and proton gyro- 
radii, the transport of heat and momentum perpendicu- 
lar to the magnetic field lines is expected to be strongly, 
but not necessarily fully, suppressed. This should lead to 
anisotropic heat and momentum transport preferentially 
along the field lines. We discuss the effect of magnetic 
fields and thermal conduction on the KHIs at the CF in 
the next subsection, here we discuss the implications of 
a possibly more complex nature of viscosity. 

In the anisotropic case, perturbations perpendicular 
to the magnetic field lines can be unstable. The insta- 
bilities are only suppressed along the field lines. In a 
tangled field geometry, however, the suppression of the 



instabilities is mostly isotropic. Dong & Stone (2009) 
demonstrated these effects for buoyant cavities in the 
ICM. Thus, the effective viscosity relevant for the slosh- 
ing CFs will dep end on the field geo metry. Simulations 
of MHD sloshing ( Zu Hone et al.|201l| ) demonstrated that 
sloshing does not fully erase the field tangling. Moreover, 
the ordered component of the magnetic fields would be 
preferentially aligned with the CFs. Given the merger ge- 
ometry in Virgo, we expect the ordered magnetic fields 
mostly in the plane of the sky. Consequently, we would 
observe mostly the KHI modes along the magnetic field 
lines, and our conclusions regarding the ability of the 
viscosity to suppress the observable KHIs should remain 
valid. 

The variations in local magnetic field structure could 
vary the local effective viscosity, and we need to consider 
to what extent the effective viscosity inferred at the CFs 
is representative f or the ICM as a whole. The simulations 
of ZuHone et al. (2011) indicate that the field amplifica- 
tion is not restricted to a narrow layer close to the fronts, 
but seems to extend throughout the core region. There- 
fore the effective viscosity inferred from the CF structure 
may be representative for the entire cluster core region. 
The relation of the viscosity inferred at the CF to the vis- 
cosity of the ICM well outside the sloshing region may 
differ, however. 
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Fig. 9. — Mock brightness residual maps for different lines-of-sight (LOS). The mock X-ray images are divided by their azimuthally 
averaged image. The left column shows the low viscosity case, the right column the high viscosity case. The top row is for the LOS 
perpendicular to the merger plane. In the subsequent rows the angle between the LOS and the normal to the orbital plane is 30, 60, and 
90 degree. The rotation is around the vertical, i.e. N-S axis. For rotation around the E-W axis see Fig. [To] 

At high viscosity, the outer CF appears as a smooth arc with a sharp edge, whereas the CF appears ragged and fuzzy at low viscosity. 
While this statement is tr ue for all LOSs, the ov erall morphology of the observed brightness residuals in Virgo agree with the top two rows 
only (compare to Fig. 1 in Roedige r et al.|2011| >. 




Fig. 10. — Same as Fig. [9] but rotation around the horizontal, i.e. E-W axis. 
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6.2. Alternative processes suppressing the KHI 

Apart from an effective viscosity, aligned magnetic 
fields or thermal conduction may lead to a suppression 
of the KHI at CFs. We discuss to what extent these 
processes can be distinguished observationally. 

6.2.1. Magnetic fields aligned with the cold fronts 

Sufficiently strong magnetic field s parallel to the shea r 
flow interface suppress the KHI ( |Chandrasekhar 1961). 
Strongly magnetized magnetic draping layers have been 
shown to form along merger cold fronts, i.e. along the 
leading edges of subclusters or galaxies moving through 
the magnetized ICM of their hos t clusters flLyu tikov 
2006; Du rsi fc Pfrommer| [20081 see |Vikhlinin & Marked 
vitch 2002] for an application to the merger CF in Abell 
3667). Such merger CFs are the contact discontinu- 
ities between the cooler atmosphere of the subcluster or 
galaxy and the hotter ICM of the host cluster. 

At first glance one might expect that a similar mag- 
netic draping layer forms at sloshing cold fronts as well, 
and could suppress the KHI. However, the mechanism 
that forms sloshing CFs is fundamentally different from 
the one responsible for the merger CFs. It is not related 
to a body moving through an ambient gas, but it is a 
nonlinear wave phenomenon of the ICM in the cluster 
potential (Nulsen & Roediger, in prep.). As a conse- 
quence, gas sloshing leads to magnetic field patterns sig- 
nificantly different from the draping layers at the merger 
CFs. Sloshing of a magnetized ICM leads to an ampli- 
fied and ordered magnetic field on the inside, the colder 
side, of the sloshing CFs, not on the outer and hotter 
side (see MHD simulations of |ZuHone et al.||201l| ). The 
amplification and "combing" or the held is not restricted 
to thin layers near the front, but layers of amplified fields 
can be found at the CFs and elsewhere on their inside. 
These layers are not necessarily thin, but can be tens of 
kpc wide. The draping layers at the merger CFs finally 
saturate at equipartition with the thermal pressure, but 
the sloshing-amplified magnetic field strength depends 
on the assumed initial field strength. While the average 
field strength throughout the cluster core has a some- 
what weaker dependence on the initial strength, the field 
strength directly at the CFs is sensitive to the initial 
strength ( | ZuHone et al. 2011| ). Consequently, the abil- 
ity of the amplified field to suppress the growth of the 
KHI depends on the initial magnetic field strength, and 
only the strongest plausible initial magnetic field reached 
amplifications sufficient to suppress the KHI, which is es- 
p ecially true for the out er CFs. 

ZuHone et al.l (|2011|) emphasized that their simula- 



netic field strengths are observationally constrained. 

Ob serv ations show both ragged and smooth CFs (see 
Sect. 6.3). In the case of smooth CFs, it would be help- 



tions demonstrated the ability of observationally plausi- 
ble magnetic fields to suppress the KHI at sloshing CFs. 
We would suggest a somewhat different interpretation of 
their simulation results: Given that only the strongest 
plausible fields clearly suppressed the KHI, which is es- 
pecially true for the outer CFs, this indicates that ob- 
served magnetic field strengths in clusters are in a range 
where they may or may not be able to suppress the KHIs 
at sloshing CFs. The answer may differ from cluster to 
cluster, and ind eed b oth smooth and ragged CFs are ob- 
served (see Sect. 6.3). A detailed investigation of individ- 
ual clusters is required, where the actual merger history, 
and thus dynamical ICM properties, as well as the mag- 



ful if we could distinguish whether the suppression of 
the KHIs occurred by viscosity or amplified magnetic 
fields. If a thin layer of magnetic fields aligned with the 
CF suppresses the KHI, the field strength in this layer 
needs to e xceed about 10% of the therm al pressure in 
this layer ( Vikhlinin fc Markevitch||2002 ) , which results 
in an approximately 20% surface brightness depression, 
which could be detected in deep observations. However, 
the sloshing-amplified fields do not have a simple draping 
layer geometry, but the amplified layers may be rather 
wide and extend far below the CFs. Consequently, the 
resulting reductions in thermal ICM pressure would not 
be localized at the CF, but the reduction would be grad- 
ual from the CF inwards, and it may be challenging to 
recognize the presence of this gradual reduction in ther- 
mal ICM pressure pr ofiles. However, the simulations of 
ZuHone et al. (2011) suggest that sloshing forms mul- 
tiple layers, streaks of amplified magnetic fields below 
the CFs. Their magnetic pressure leads to corresponding 
streaks of reduced ICM density that should be detectable 
as ripples in surface brightness. It will be interesting to 
compare not only the structure at the CFs, but also below 
the CFs in order to distinguish the cases of high viscos- 
ity and amplified aligned magnetic fields. Additionally, 
a significantly amplified field over the whole cluster core 
may well give rise to radio synch rotron emission and thus 
radio minihalos, as s uggested b y Mazz oTta fc Giacintucci| 
fl2008|) and shown by |ZuHone et al.l Q201ip . It will belm 
teresting to investigate the correlation between the pres- 
ence of radio minihalos and the smoothness of sloshing 
CFs. 

6.2.2. Magnetic fields and anisotropic viscosity 

It will be interesting to study the combined effect of 
magnetic fields and (anisotropic) viscosity, because they 
have opposite effects and mutual dependencies. In the 
standard picture, zero magnetic fields imply a high vis- 
cosity, and thus sharp CFs that form smooth arcs. In- 
creasing the magnetic field strength is expected to reduce 
the effective viscosity on the one hand, thus permitting 
KHIs larger than a certain wavelength. Thus, as a first ef- 
fect, the presence of magnetic fields facilitates the growth 
of KHIs. This could be observed as a ragged structure 
of the CF arc, and the density and brightness jump be- 
ing washed out. The instabilities create turbulence and 
may tangle and amplify the fields, reducing the effec- 
tive viscosity further. However, once the magnetic fields 
are strong enough, they may, themselves, suppress the 
growth of the KHI. Their effect on the ICM density and 
temperature structure should be very similar to the vis- 
cous suppression of the KHI, but the enhanced magnetic 
fields may be revealed by the radio minihalos. 

6.2.3. Thermal conduction 



Xiang et al. (2007) discuss the impact of thermal con- 
duction on merger CFs. Given that the formation mecha- 
nism of merger and sloshing CFs differs, the dire ct appli- 
cation of their results to sloshing CFs is unclear. ZuHone 
et al. ( |2012[ ) presented MHD simulations of gas sloshing 
and included an anisotropic thermal conduction, where 
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conduction was fully or strongly suppressed perpendic- 
ular to the magnetic field lines. The sloshing leads to 
an approximate alignment of the magnetic fields with 
the fronts, but the cool regions below the fronts are still 
connected with hotter regions above the fronts, and full 
Spitzer conduction along the field lines has a significant 
impact on the CFs. The temperature and density con- 
trasts are strongly reduced, and the fronts are less sharp, 
the discontinuity is washed out. The resulting finite 
width of the interface suppresses KHIs of perturbation 
length scales about several t imes as long as the interface 
width ( |Chandrasekhar|1961[ ) . 

In contrast to viscosity, thermal conduction widens the 
temperature and density discontinuity, whereas viscosity 
widens the discontinuity in shear velocity. This leads 
to a characteristic difference in the CF structure: Ther- 
mal conduction should lead to CFs free of KHI below 
a certain wavelength, but a washed out discontinuity in 
density and temperature, i.e. the edge will not be sharp 
in either of these two quantities. Viscosity, on the other 
hand, will lead to a smooth-arc CF as well, but a sharp 
one, as neither the density nor the temperature disconti- 
nuity are washed out. The presence of KHIs leads to an 
apparently washed-out CF in surface brightness as well, 
but in this case, individual KH rolls may be identifiable. 
Moreover, the cold gas dominates the emission measure 
in the mixing layer induced by the KHIs, and we expect a 
steep jump in projected temperature across the front. In 
contrast, thermal conduction truly washes out the tem- 
perature jump, leading to a different profile in projected 
temperature across the front. 

If either thermal conduction or viscosity suppresses the 
KHI not fully but only below an observable length scale, 
KH rolls of larger length should be recognizable. As just 
mentioned, if viscosity suppresses the small KHIs, in- 
dividual stretches of the CF shorter than the limiting 
length scale should appear as sharp edges, and the ex- 
isting, larger KH rolls may be very pronounced (see the 
big KH rolls in the SW of the simulated Virgo cluster 
in Fig. [3|. In the case thermal conduct ion suppresses 
the sm all-scale KHIs, the simulations of ZuH one et al.| 
( |2Q12| ) indicate that short stretches of the CF should not 
be sharp but be washed out. 

The CF structure resulting from the combined effect 
of thermal conduction and viscosity has yet to be deter- 
mined. 

6.3. Comparison to other clusters 

The detailed structure of sloshing CFs differs be- 
tween clusters. In some clusters they appear as smooth 
arcs within the observational resolution, e.g. Abell 2142 
(iMarkevitch et al | | 2000] |Owers et aL]|2QQ9| ) and Abell 
2204 (ganders et al.||2UU9jT The interpr etation of the 
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in Abell 496 (Dupke et al. |2007| ) have a distinct boxy or 
even concave app earance tnat resembles KH rolls ( |Roedi-| 
ger et al. 2012b). The CF towards the west is a dou- 
blet similar to the multiple adjacent edges reported here. 
The northern CF seems to curve around in the NE, but 
becomes washed out in the east. The sloshing CFs in 
the group-group merger NGC 7618/UGC 12491 display 
a triangular-shaped nose and wings reminiscent of KH 



rolls shown in Fig. [3] ( |Roediger et al. 2012a). We note 
that a significant fraction oi the observed cold fronts in 
cluster cores are not shaped like the long arcs in spiral 
or staggered arrangements predicted by numerical simu- 
lations, but several clusters exhibit only short stretches 
of clear cold fronts. It will be interesting to revisit the 
observational data and attempt to distinguish whether 
this is due to the orientation of the merger with respect 
to our LOS, an intrinsic ellipticity of the cluster or a 
washing out of the CFs by instabilities. 

7. SUMMARY 

In this paper we present the first detailed attempt to 
constrain the effective viscosity of the ICM from the pres- 
ence or absence of KHIs at sloshing CFs. We focussed on 
the Virgo cluster for its proximity, brightness and favor- 
able merger geometry. We performed viscous hydrody- 
namical simulations tailored for the Virgo cluster, com- 
pared them to existing XMM-Newton observations and 
made predictions for future deep Chandra observations. 
Our key results are: 

• At the Virgo cluster cold fronts, a viscosity of 10% 
of the Spitzer value suppresses the KHI for scales 
below tens of kpc. In contrast, the KHI occurs 
at ~ 15 kpc scales at viscosities below 1% of the 
Spitzer value (Fig. [2]). 

• We identify the outer CF at ~ 90 kpc to the north 
and north-east of the Virgo center as the best tar- 
get for distinguishing different viscosities observa- 
tionally, because this front is not disturbed by nu- 
clear outbursts of M87, and it is subject to a suffi- 
cient shear flow which potentially induces the KHI. 

• Our simulations predict observationally distin- 
guishable structures at this CF: at viscosities of 
> 10% of the Spitzer value, the KHI is suppressed 
along the northern CF.it appears as a smooth arc 
(bottom panel of Fig. pi). Surface brightness pro- 
files across the front are consistent with a sharp 
density discontinuity (bottom panel of Fig. [7|. At 
viscosities < 1% of the Spitzer value, the front has 
a ragged, sawtooth-like appearance (top panel of 
Fig. [3). Individual KH rolls at ~ 15 kpc size can be 
identified, and the front splits up into multiple, ad- 
jacent brightness edges. Surface brightness profiles 
across the front are shallower than expected for a 
density discontinuity, i.e. the discontinuity appears 
blurred out. Profiles in narrow wedges show mul- 
tiple jumps arising from individual KH rolls (top 
panel of Fig. [7]). 

• We re- analyze archival XMM-Newton data for cor- 
responding signatures. The northern cold front ap- 
pears to be less sharp in the north-east than in the 
north-west (Fig. [5]), indicating a blurring out of the 
discontinuity in the north-east by KHIs and thus 
a low ICM viscosity. Due to low count statistics 
and the non-uniform detector response these data 
do not allow the detection of individual KH rolls. 

• A 300 ks Chandra observation of the north-eastern 
front would confirm the presence or absence of 
KHIs both in images (Fig. m and in profiles across 
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the front (Fig. [7| and thus constrain the effective 
ICM viscosity. 

• We discuss in detail the differences in KHI sup- 
pression by viscosity, magnetic fields and thermal 
conduction, and how these processes can be distin- 
guished observationally. 
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APPENDIX 



KELVIN-HELMHOLTZ INSTABILITY AT LOW RESOLUTION 
Figure [TT] demonstrates the ability of the FLASH code to capture the KHI even at low resolution. 
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Fig. 11. — Resolution test for the viscous Kelvin-Helmholtz instability in an idealized setup. The simulation box is periodic in x-direction. 
Pressure is constant throughout the simulation box. Above x = the fluid temperature is twice as hot as below. The upper fluid moves 
right, the lower fluid moves left, their relative velocity is Mach 0.5. Initially the interface is perturbed with a sinusodial vertical velocity v y 
of amplitude Mach 0.05 and a wave length of half the simulation box size. We assume a Spitzer-like viscosity with a suppression factor set 
such that the Reynolds number in the hotter layer is 100. See Roediger et al., in prep., for more details. 

The panels show the density of a tracer fluid that was set to 1 in the lower layer and zero in the upper. Rows from top to bottom are for 
resolutions of 128, 64, 32, 16 cells per perturbation wavelength. Left, middle and right columns are for timesteps of 50 Myr, 100 Myr, 200 
Myr (6.25, 12.5, 24 growth times). The FLASH code captures the instability even at the lowest resolution. 



